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Abstract 

Gribov approach to high-energy interactions of hadrons and nuclei is reviewed 
and applied to calculation of particle production in heavy-ions collisions. It is 
pointed out that the AGK (Abramovsky, Gribov, Kancheli) cutting rules is a po- 
werful tool to investigate particle spectra in these processes. It leads, in the Glauber 
approximation, to a simple formula for the density of hadrons produced in the cen- 
tral rapidity region in nucleus-nucleus interactions. An estimate of this density 
for RHIC and LHC is presented and compared with results of Monte-Carlo cal- 
culations. It is shown that the Glauber approximation substantially overestimate 
particle densities compared to the results of the complete Gribov theory. This is due 
to extra shadowing in the system, related to large mass diffraction which leads to 
a strong decrease of particle densities at mid rapidities. Our method of calculation 
of these effects has been applied to the problem of shadowing of nuclear structure 
functions and a good agreement with experimental data has been obtained. 
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1 Introduction 



The reggeon approach to high-energy interactions is undoubtly an important ingredient of 
modern theory. V.N. Gribov has made very essential contributions to the development of 
this approach. He has introduced a leading pole in the complex angular momentum plane 
JT|], which determines asymptotics of diffractive processes (nowadays called the Pomeron), 
investigated the main properties of Regge-poles and Regge-cuts || and developed the 
reggeon diagram technique ||. 

An important contribution to the theory of multiparticle production has been made 
by V.N. Gribov together with V.A. Abramovsky and O.V. Kancheli || and is usually 
referred to as the AGK-cutting rules. These rules are widely used for the construction 
of QCD-based models for high-energy interactions (for reviews see |6|, 0]). In this paper 
we will extensively use AGK-cutting rules and Gribov reggeon diagram technique to 
calculate inclusive particle spectra and particle densities in the central rapidity region 
for heavy-ions collisions at future colliders RHIC and LHC. The values of particle and 
energy densities are very important for the problem of creation of the quark-gluon plasma 
in heavy-ions collisions at these very high energies. 

In the first part of the paper (Section 2) we will shortly review Gribov's approach to 
high-energy interactions of hadrons with nuclei || and note the difference in the space- 
time picture of interactions at low and high energies. We will also discuss a difference be- 
tween the Glauber approximation and the general Gribov theory. The AGK-cutting rules 
will be applied to hadron-nucleus and nucleus-nucleus collisions and a simple formula for 
inclusive particle spectra and densities of produced particles in the central rapidity region 
valid in the Glauber approximation will be presented. 

Applications of the Regge-Gribov theory to high-energy hadronic interactions are 
considered in Section 3. The important role of shadowing effects in these processes is 
emphasized. It is pointed out that understanding diffractive processes is important for 
a selfconsistent description of high-energy interactions of hadrons and nuclei. The re- 
sults of this analysis are used in the Section 4 to estimate particle densities at RHIC 
and LHC . Limitations of the Glauber approximation and its modification at high energy 
according to the Gribov theory will be discussed . It will be shown that these modifi- 
cations are related to large mass diffraction dissociation of hadrons which leads to extra 
shadowing in the system. These effects reduce particle densities in the central rapidity 
region compared to the Glauber approximation predictions. This result is valid for both 
soft and hard processes. To estimate these shadowing effects we apply Gribov theory to 
the processes of deep inelastic scattering (DIS) on nuclei and show how it is possible to 
calculate nuclear shadowing effects in the region of small Bjorken-x, using information on 
diffractive production in DIS. Our results agree with existing experimental data and allow 
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a safe extrapolation of shadowing effects to the experimentally unmeasured region of very 
small x, relevant for nuclear collisions at LHC This model leads to a reduction of particles 
densities in PbPb collisions at LHC (RHIC) by a factor of three (two), with respect to 
the results of the Glauber approximation. Our results are compared to predictions of 
different Monte-Carlo calculations. 



2 Gribov theory of high-energy nuclear interactions 

In classical papers of Gribov || it was shown how to incorporate the Glauber model || for 
interactions of hadrons with nuclei into a general framework of relativistic quantum theory. 
Consider an amplitude of elastic scattering for high-energy hadron-nucleus interactions. 
In the Glauber model it is described by the diagrams shown in Fig. 1, which looks like 
a successive rescatterings of initial hadron on nucleons of the nucleus. However, as was 



emphasized by Gribov || |Kj the space-time picture of the interaction at high-energy 
E > m^nR-A (/•* is a characteristic hadronic scale ~ 1 GeV and Ra is the radius of 
the nucleus) is completely different from this simple picture. It corresponds to coherent 
interactions of a fluctuation of the initial hadron, which is "prepared" long before its 
interaction with the nucleus (Fig. 2). A very important result of Gribov || is that 
nevertheless the elastic hA-amplitude can be written as a sum of the diagrams shown in 
Fig. 3, with elastic rescatterings (Fig. 3a) which give the same result as Glauber model, 
plus all possible diffract ive excitations of the initial hadron. At not too high energies 
El ~ 10 2 GeV these terms lead to corrections to the Glauber approximation of 10 — 20% 



for the total hA cross section [11, 12 



We will show below that at very high energies and for inclusive cross sections this 
modification of the Glauber approximation is very important. The difference between 
Glauber model and Gribov's theory is essential for understanding shadowing corrections 
for structure functions of nuclei related to hard processes on nuclei and for many aspects 



of multiparticle production on nuclei [JT3 



An important consequence of the space-time structure of the diagrams of Fig. 2 for 
interactions of hadrons with nuclei is the theorem, based on AGK-cutting rules ||, that 
for inclusive cross sections all rescatterings cancel with each other and these cross sections 
are determined by the diagrams shown in Fig. 4 (impulse approximation). Note, however, 
that this result, valid asymptotically in the central rapidity region, only applies to the 
diagrams of the Glauber-type, i.e. when masses of intermediate states in Fig. 3 are limited 
and do not increase with energy. As a result, the inclusive cross section for the production 
of a hadron a is expressed, for a given impact parameter b, in terms of inclusive cross 
section for hN interactions 
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d 3 p ' d 3 p 

where T A {b) is the nuclear profile function (/ d 2 b T A (b) = A). After integration over b we 
get 

^ d 3 af lA _ ^ g d 3 <JhN 
d 3 p d 3 p 

The total and inelastic hA cross sections in the Glauber model can be easily calculated 
and are given for heavy nuclei by well known expressions. For example 

at, = J d 2 b(l - exp(-a™ N T A (b)) (3) 

The situation for nucleus-nucleus collisions is much more complicated. There are no 
analytic expressions in the Glauber model for heavy-nuclei elastic scattering amplitudes. 
The problem stems from a complicated combinatorics and from the existence of dynamical 
correlations related to "loop diagrams" |EI| [HJ . Thus, usually optical-type approximation 
16, [L7| and probabilistic models for multiple rescatterings fl8| are used. For inclusive 



cross sections in AB-collisions the result of the Glauber approximation is very simple to 
formulate due to the AGK cancellation theorem. It is possible to prove, for an arbitrary 
number of interactions of nucleons of both nuclei fll9| , that all rescatterings cancel in 



the same way as for hA-interactions and only the diagrams of Fig. 5 contribute to the 
single inclusive spectrum. Thus a natural generalization of eq. ([3]) for inclusive spectra 
of hadrons produced in the central rapidity region in nucleus-nucleus interactions takes 
place in the Glauber approximation 



E U U ff KU) = T AB {b) (4) 



,d 3 cr a AB (b) _ E fa a 
d 3 p ' 1 ' d 3 p 

where T A B{b) = J d 2 sT A (s)TB(b — s) . After integration over b eq. (§) reads 



E d^B = AB E d L Oj m _ 
d A p d A p 

Note that eqs.([|),(]5]) are valid for an arbitrary set of Glauber diagrams and are not 
influenced by the problem of summation of "loop" diagrams mentioned above. 

The densities of charged particles can be obtained from eqs.([|), (|5|) by deviding them 
by the total inelastic cross section of nucleus-nucleus interaction. For example 

dn%{b) = T AB (b)da% h N 
dy (?ab dy 
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and 

dn c AB _ AB da c NN 
dy ~ a% dy 1 ' 

In the following we shall use these results to calculate particle densities in the central 
rapidity region at energies of RHIC and LHC. 

3 Regge-Gribov theory and hadronic interactions. 

In this section we will briefly review results obtained as an application of the Gribov 
theory to hadronic interactions at very high energies. We will need them in order to 
determine inclusive particle spectra in NN-interactions at RHIC and LHC and in order 
to illustrate the importance of inelastic diffractive processes for a selfconsistent treatment 
of high-energy hadronic interactions. 

In the traditional Gribov's reggeon approach it is assumed that the Pomeron is a 
Regge pole accompanied by the cuts associated to the multi-Pomeron exchanges in the 
t-channel. These cuts are important to restore the unitarity of the theory. Reggeon 
diagrammatic technique and the AGK cutting rules allow one to calculate con- 
tributions of many-Pomeron exchanges to scattering amplitudes and relate them to the 
properties of multiparticle production. 

An important parameter of the theory is the value of the Pomeron intercept ap(0). 
In the simplest model, where only the single-Pomeron exchange is taken into account the 
intercept determined from the analysis of a tot and elastic scattering data is ap(0) ~ 1.08 
20[| . However, this is only an effective value of the intercept a e J (0), which describes 



the energy dependence of total hadronic cross sections in the currently available range 
of energies, - er^ *) ~ s Q p // (°)- 1 , An extensive phenomenological analysis, which takes 
into account multi-Pomeron exchanges, (see e.g. refs. [[], 0]) shows that the Pomeron 
intercept is substantially larger than the value ap (0) indicated above. With eikonal type 
diagrams only one gets ap(0) = 1.12 -r- 1.15 For a more complete set of diagrams, 
which include interactions between exchanged Pomerons (related to large mass diffractive 
production), an even larger intercept of ap(0) ~ 1.2 is obtained f22f . 

In this approach many characteristics of high energy hadronic interactions are well 
described |6|, [7J. Multi-Pomeron exchanges are very important for understanding many 
qualitative features of experimental data. For example the fast increase of inclusive spectra 
in the central rapidity region can be reproduced only if multi-chain configurations, which 
are due to cutting of multi-Pomeron diagrams, are taken into account. For pure pole 
model the density of produced hadrons in the central region would be energy independent 
at large s. Experimentally in pp-interactions it increases with energy approximately as 
s s with 5 « 0.13. 
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Another manifestation of multipomeron exchanges is an existence of important long- 
range correlations (for example forward-backward correlations). They are firmly estab- 
lished experimentally at high energies. Long range correlations are closely related to 
broad multiplicity distributions. The models based on the reggeon diagrams technique 
and AGK-cutting rules give a good quantitative descriptions of multiplicity distributions 
including violation of the KNO-scaling at very high energies. 

Let us consider now diffractive production of hadrons at very high energies. Descrip- 



tion of these processes in the Regge-model can be found in reviews p3| , |24| , |25| . The 
differential cross section for inclusive single diffraction dissociation can be written in the 
form [g3j 

d 2(y (9pp{t)) 2 , tl ,x,2 {tot), ,s / Q x 

where S2 = M 2 , M is an invariant mass of the produced state, Gp(£',t) = 
T](ap(t)) exp[(ap(t) — 1)£], £ = ln(s/s2) and rj(ap(t)) is the signature factor. The quantity 
a f pp(s2,t) is the total cross section for Pomeron-particle interaction. This cross section 
has asymptotic Regge behavior for large mass values of s 2 - the squared of diffractively 
produced state. In this case, the cross section for diffraction dissociation is described by 
the triple-Regge diagrams (Fig. 6) and has the form 



^^i^, t ),*p>)*(r (0 " 

Values of <Tp^\s^ t) and of the triple reggeon vertices rp P , r pp have been determined 



from analyses of diffractive production in hadronic collisions || ^4], ^ . 

In the pole approximation with A = 0.08 the total cross section of inelastic diffrac- 
tion increases too fast with energy and strongly contradicts to experimental data at 
y^s ~ 10 3 GeV. This problem is solved by the inclusion of Regge cuts. For example 
in ref. a much slower increase of <td, consistent with recent experimental result, was 
predicted even for A = 0.21 . It is important to note that in Gribov theory the amount 
of shadowing is closely related to the magnitude of diffractive processes, which in their 
turn are influenced by the shadowing. This complicated, nonlinear problem can be solved 
only by a systematic and selfconsistent treatment of both diffractive processes (elastic 
and inelastic) and multiparticle production. 
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4 Particle densities in heavy— ions collisions at super- 
high energies 

Now we will address the question of particle densities in heavy-ions collisions at energies 
of future colliders. For central PbPb collisions at LHC (y/s ~ 6 TeV), existing Monte- 



Carlo models (see [27]) give a rather broad range of values, of dn jdy at y* ~ -from 



about 8000 particles for the VENUS model [28] to ~ 1400 particles for the String Fusion 



Model (SFM) |29|| . Other Monte-Carlo models j3(| give predictions within this inter- 
val. All these models are based on the probabilistic approximation to the Glauber model 
of nuclear interactions and should lead to similar results for inclusive particle densities, 
which according to eqs.(|6]) -([?[), do not depend on details of the particular model (for NN- 
interactions all these models give similar predictions, since the extrapolation from the 
experimentally measured region of energies is not large). Below, we shall compare these 



results of Monte-Carlo models and of semi-analytic calculations |32| with predictions of 
eqs.(|),(^) and shall discuss deviations from the Glauber approximation build in some 
of these models. Finally, we will compare them with predictions of the complete Gri- 
bov theory of nuclear collisions -i.e. including shadowing corrections due to high mass 
intermediate states in Fig. 3. 

Eq. (0) for particle densities integrated over impact parameter (minimum bias events) 
can be rewritten as 

dn c AB dn c ^ N 

~dy~ = B ~dy~ (10) 

where n^B — — sr^- It corresponds to the average number of collisions in the Glauber 

model. For A = B » 1 n AB behaves as CM 4/3 with C w ^§ (R A = RoA 1 / 3 ). It 
is well known that eqs.(P|),(|5|),(|lCi|) can be applied to hard processes but in the Glauber 
approximation they are valid for soft processes as well. We shall see below that for both 
soft and hard processes these equations have to be modified. 

For characteristics of pp-interaction at y/s = 6 TeV we take predictions of the Quark- 
Gluon String model and dual parton model (DPM) |6|, |7j : total inelastic cross section 

= 65 mb, inelastic nondiffractive cross section a l Jl ,nonD = 50 mb, and for density of 
charged particles at y = for nondiffractive interactions drf^jdy = 5.0. The uncertainty 
in these numbers is ~ 10%. The total inelastic cross section for nucleus-nucleus collisions 
can be calculated either using a simple geometrical formula or the optical approximation 
to Glauber model, which both give at LHC energies for PbPb collisions a value ~ 5 barn 
(in these estimates an increase of the radius of NN-interaction with energy has been 
taken into account). Using these numbers and eqs. (|7|), (|ToD we obtain for PbPb collisions 
at LHC at y* = the following numbers for minimum bias events and central (b < 3 fm) 
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collisions, respectively 

dn ch /dy = 2100 , dn ch /dy = 8500 . (11) 



These numbers are close to results of the VENUS model pSJ and to the DPM results 
of ref. [B2] (1890 and 7900 respectively, with a slightly different definition of central 



collisions). However, the results of SFM |29| without string fusion and DPMJET-II [31 



are about twice smaller. For SFM this can be due to the fact that this code has limits on 
the number of produced strings and is not reliable at LHC energies JEJ . 

Thus, the Glauber approximation predicts very large densities of charged hadrons in 
central heavy-ions collisions at LHC. However, are these predictions realistic ? In order 
to answer this question we will consider possible limitations of the Glauber approximation 
and also the corrections to the AGK cancellation theorem which are important at high 
energies. 

There are two types of corrections to eqs.(]I|),(|j), 

a) The effects due to energy-momentum conservation [34|, -the energy of the initial 
hadron is shared by "constituents" (see Fig. 2) and each sub-collision happens at smaller 
energy. These effects are very important in the fragmentation regions of colliding hadrons 
(or nuclei) and reduce particle densities. For y* = this reduction decreases as (l/s) 1 / 4 . 
It is important at SPS energies ; however, at LHC energies in the central rapidity region 
this effect is small. It is taken into account in the Monte-Carlo models mentioned above. 

b) Another dynamical effect is important at very high energies when diffractive production 
of very heavy hadronic states (M 2 >> m 2 N ) becomes possible and should be taken into 
account in the diagrams of Fig. 3. Consider for example a double rescattering diagram of 
a proton on a nucleus, which contains, according to Gribov theory, the diffractive large- 
mass intermediate states shown in Fig. 7. It is related to the triple-Pomeron interaction 
discussed above and corresponds to an interaction between Pomerons (strings in the string 
models of particle production). As the total and inelastic cross sections of hA and AB- 
interactions at high energies are close to a black disc limit due to Glauber-type diagrams, 
these extra interactions have a small influence on total cross sections. However, they are 



very important for inclusive spectra in the central rapidity region ||13|| , where contributions 
of Glauber rescatterings cancel due to AGK rules. 

The diagram shown in Fig. 7 is only one of a large class of diagrams with interactions 
between Pomerons. The application of AGK-cutting rules to these diagrams leads to the 
diagrams for inclusive cross section in AB-collisions shown in Fig. 8. Extra shadowing 
effects related to these diagrams modify the A-dependence of the Glauber approximation 
for inclusive spectra eqs.(|l]),([D in such a way that the behaviour ddAB/dy ~ AB of the 
Glauber approximation changes to da^/dy ~ A a B a , where a < 1. For very strong 
interaction between Pomerons a — > 2/3 . This limit was considered by O.V. Kancheli a 
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long time ago It leads to universal particle densities in pp, pA and AB-collisions. We 
will show that due to a rather week interaction between Pomerons even at LHC energies 
the value of a is close to 0.9. 

The problem of shadowing for inclusive spectra is not especially related to soft pro- 
cesses. The same graphs of Fig. 8 are relevant also for hard processes (production of 
jets or particles with large pr, heavy quarks, large-mass lepton pairs, e t.c). For hard 
processes, due to the QCD-factorization theorem, inclusive spectra in nucleus-nucleus 
collisions are given by convolutions of hard cross sections with distributions of partons in 
the colliding nuclei. In these cases the diagrams of Fig. 8 describe shadowing effects for 
nuclear structure functions (i.e. distributions of quarks and gluons in nuclei). Due to a 
coherence condition these effects are important only in the region of very small x,- t of par- 
tons ( Xi « 1/ {Rattin))- So these effects are important only at very high energies, when 
Xi ~ Mj-j^fs satisfy this condition. This condition in terms of Xj of partons coincides 
with the condition on diffractive production of large-mass states discussed above (see for 
example fl36| ). 



The effects of shadowing are observed experimentally in deep inelastic scattering on 



nuclei |37], |38|. So in order to test the Gribov method of calculations of shadowing 
corrections we will apply it to these processes (Fig. 9). The first diagram in Fig. 9 
corresponds to a sum of interactions with nucleons of nuclei and is proportional to A. The 
second diagram (Fig. 9b)) describes the shadowing effect due to a coherent interaction of 
a virtual photon with two nucleons of the nucleus and is related to diffractive production 
in DIS. This process was measured at HERA |3£], [fij and was well described in the model, 
based on the triple-Regge diagrams of Fig. 10 |jfl, |42|. The study of diffraction dissociation 
of a virtual photon at HERA allows a better determination of triple-Regge couplings 
compared to hadronic reactions. This is related to the fact that diffraction dissociation 
in DIS is much less influenced by absorptive corrections than diffractive production in 
hadronic interactions, where "effective" vertices are much smaller than their "bare" values 
(see below, after eq. (|T7|)). 

The contribution of a double rescattering term to the directly expressed in 

terms of the differential cross section for the diffraction dissociation of a virtual photon 
in 7*AMnteractions: 

^ (2) = / d%Tl(b) J dM* da ^l = 0) ] F A (t mm ) (12) 

where FA{t m in) = exp(.R^t m j n /3); t m i n ~ —m? N x 2 p. 

Higher order rescatterings are model dependent and in the generalized Schwimmer 
model we obtain for the ratio F^aIF^n structure functions F 2 nuclei and nucleons, in 
the region of small x 
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^W ^ + ^IVw (13) 

with F(i,Q 2 ) = 4ir/dAf 2 (<fo£ft(t = 0)/dM 2 (i() (F^(i„i„)/<T 7 -jv(a:,<3 2 )). 

Theoretical predictions |0J for nuclear shadowing based on the model for diffraction 



dissociation of ref. |42| are in a good agreement with NMC-data on nuclear structure 



functions at very small x |37[| , including very accurate data on Sn target (Fig. 11). The 
model describes both the x and Q 2 -dependence of structure functions of nuclei and allows 
us to obtain reliable predictions for structure fmctions of nuclei in the region of x < 
10~ 3 (Fig. 12) relevant for nuclear interactions at LHC energies. Let us mention at this 
point that while the distributions of quarks in nuclei are known experimentally and can 
be rather safely calculated in the experimentally unmeasured region of x and Q 2 , the 
situation with gluons is much less clear. This problem is related to the distribution of 
gluons in the Pomeron, which can be in principle extracted from experimental data on 
diffractive production of jets or heavy-quarks in DIS. However, existing experimental data 
do not allow a reliable determination of this contribution. The terms with triple-Pomeron 
interaction, shown in Figs. 8,10 lead to a universal shadowing. So in the following we will 
assume that the shadowing effects at very small x are the same for quarks and gluons (the 



same assumptions was made in refs. |45j while in some papers |S| |47j it is assumed 



that, due to a larger ggg-coupling compared to ggg-coupling, the shadowing for gluons is 
larger than for quarks). 

The sum of diagrams of Fig. 5 and Fig. 8 leads to the following expression for inclusive 
spectra in nucleus-nucleus interactions. 

E^(b) = J dhf A (b)f B (s- b)E^L (14) 

where, in the Shwimmer model, the function /a(^) coincides with i + p^x A (b) introduced 
above in the calculation of shadowing for nuclear structure functions. After integration 



over impact parameter b eq. ( |I4"D becomes 



= F A (s A )F B (s B )E^ (15) 

Where function F A is the same as F2a/F 2 n defined above and is equal to A when the triple- 
Pomeron interactions is switched off. At y* = 0, sa = s B = rn T \J~s (which corresponds to 
Xi = 774/ y^) . 

Thus, due to interactions between Pomerons, the Glauber formula ( |Tl"| ) is modified in 
a simple way 

dn C A B dn AB Glaub. 

Hy~ = -^ly— lAlB (16) 
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where the quantity — AB is the Glauber approximation result -(eq. (0)) and 7,4 is 
the shadowing correction to nuclear structure function F 2 aI ' AF 2 n shown in Fig. 12. In 
the following, we will consider it at a low scale Q 2 ~ 1 GeV 2 relevant to our problem. For 
LHC energies 7p& ~ 0.5 — 0.6 and the total correction to the Glauber approximation for 
PbPb-collisions is 0.25 -j- 0.36. Formula ([H]) is written for minimum bias collisions. The 
result for different values of impact parameter b can be calculated using eq. fll4]) . The 
calculation shows that for central collisions this correction differs by less than 10% from 
the value indicated above. Thus there is a suppression by a factor of ~ 3 for charged 
particle densities at LHC energies compared to the results of the Glauber approximation 
given in (|ll|). For PbPb-collisions at RHIC, particle densities are reduced by a factor 
w 2 compared to the DPM predictions P2| , which do not include shadowing corrections. 



However, in this case there are large uncertainties in the calculations of these corrections 
due to a strong dependence of the quantities ^a on x« in the region of %i ~ 10~ 2 , relevant 
for these energies. 

Let us emphasize that although our derivation of the shadowing corrections is based 
on the study of nuclear structure functions, which is valid for hard processes, we do 
not assume that hard processes dominate particle production at LHC-energies. In fact, 
the shadowing corrections due to triple-Pomeron interactions are universal and the same 
correction factors in eq. ( |16D are also obtained for soft interactions by a direct calculation 
of the diagrams of Fig. 8. In the Shwimmer model we will obtain the same expression 
i+F^y)T A (b) f° r shadowing where the function F(s, y) can again be written as an integral 
of the ratio of the triple Pomeron cross-section, eq. @, over the single Pomeron exchange 
cross-section <Jp(y) : 



/Umax 
dy 



y max _ 1 d 2 a ppp 
t 

min 

« P (0)) 2 e 



a P (y) dy dt 



1 d£exp(A£) 



[exp(Ay max ) - exp(Ay min )} (17) 

where y ma x = Uc.m. + |^ n ( s / ' m2 v) > Vmin — d-n{RAm^). We use the same parameters as in 
ref. [|2| : (^(O)) 2 = 23 mb is the Pomeron-proton coupling, A = «p(0) — 1 = 0.13, e = 
rppp(0)/g%,(0) « 0.07. Eq. QTT|) leads then to practically the same suppression factors as 
given above. Note that we are using here the same value of the triple Pomeron coupling 
as in DIS calculations |42jl . As explained there, this value is about three times larger 



than the one obtained from a fit of soft hadronic diffraction using only triple Regge terms 



(without eikonal unitarization) as in refs. || [24], |25|]. Such a large value of the triple 



Pomeron coupling is required in order to describe diffraction in DIS. In ref. [42|, this was 



justified from the work of ref. ||26|| , where it has been shown that eikonalization of (Tsd 
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reduces its value by a factor close to 3 - implying that the true (or bare) triple Pomeron 
coupling is about 3 times larger than the effective one. In DIS, the eikonal corrections 
disappear very fast when Q increases and thus the bare coupling is the relevant one. In 
single particle inclusive production at high energies, the same is true since the eikonal 
corrections are absent due to the AGK cancellation. This is the physical reason for the 
large shadowing corrections obtained in this paper. 

For collisions of identical nuclei (SS, PbPb) the A 4//3 -dependence of particle densities 
of eq. (10) typical for the Glauber model changes to the behaviour A s . The value of 
delta is a weak function of energy and it is equal to 5 ~ 1.1 at LHC energies. It means 
that at these energies the values of a in the A-dependence of inclusive cross sections for 
AB-collisions daAB/dy ~ A a B a is close to 0.88. The value of a should slowly decrease as 
energy increases. In the case of stronger shadowing for gluons than for quarks a somewhat 
smaller value of a can be obtained. 

Let us compare our results with Monte-Carlo calculations, which take into account 
shadowing effects [^, [30|]. In the SFM model [2J| the interaction between Pomerons 
is introduced via a mechanism of string fusion and is estimated from geometrical sizes 
of strings. The accuracy of such an estimate is not clear but it leads to a reasonable 
suppression factor (about 2) for particle densities at LHC energies (though application 
of the SFM Monte-Carlo is questionable at LHC ,-see above). In the Hijing model |50| , 
existing data on nuclear shadowing were parameterized and thus the shadowing effects are 
also not very diffrent from our predictions, though the Hijing model leads to a somewhat 
smaller suppression. This is connected to the choice of a saturation for shadowing at 
small values of x (x ~ 10~ 3 ) in their parameterization. The reggeon formalism allows us 
to determine the x-dependence of shadowing at small x and it shows that the shadowing 
is still increasing as x decreases even for x ~ 10 -4 . Saturation happens at much smaller 
values of x. 



5 Conclusions 

Gribov theory of high-energy interactions of hadrons and nuclei is based on general pro- 
perties of amplitudes in relativistic quantum theory and provides an unified approach 
to a broad class of processes. In this theory, the Glauber approximation to nuclear 
dynamics is valid in the region of not too high energies and should be modified at energies 
of RHIC and LHC. AGK- cutting rules provide a very powerful tool for the study of 
multiparticle production for all types of high-energy processes and allow one to obtain 
simple predictions for inclusive cross sections in hh, hA and AB-collisions. 

In this paper we used AGK-cutting rules to obtain predictions for densities of particles 
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at future heavy-ions colliders. Gribov theory then allows to determine corrections to the 
Glauber approximation for inclusive particle spectra by relating them to cross sections 
of large-mass diffraction. The technique has been applied to calculation of shadowing 
effects for structure functions of nuclei and a good agreement with experimental data 
on these processes has been obtained. The same approach predicts a strong reduction 
of particle densities at super-high energies as compared to predictions of the Glauber 



approximation. Our calculations show that the DPM results for PbPb collisions are 
reduced by a factor 3 at LHC and 2 at RHIC energies. The expected values of the charged 
particle density in central PbPb collisions at y* = is thus 2500 at LHC and 1000 at 
RHIC. 

Future experiments at RHIC and LHC will test these theoretical predictions and will 
allow a better determination of the parameters that govern the dynamics of shadowing. 
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FIGURE CAPTIONS 
Fig. 1 Diagrams of the Glauber model for the elastic /M-scattering amplitude. 

Fig. 2 Diagrams for high-energy /M-interactions. 

Fig. 3 The dispersion representation diagrams for /iA-elastic scattering amplitude. A 
cross on a line means that the particle is on the mass shell. 

Fig. 4 The diagram for inclusive cross section of particle a in /iA-collisions. 

Fig. 5 The diagram for inclusive cross section in the Glauber approximation for AB- 
collisions. 

Fig. 6 Triple-regge diagrams for diffractive production of large-mass states in pp- 
collisions. 

Fig. 7 Large-mass diffractive contribution to pA-elastic scattering amplitude. 

Fig. 8 Diagrams for inclusive cross sections in AB-collisions, which take into account 
interactions between Pomerons. 

Fig. 9 Diagrams for interactions. 

Fig. 10 Triple-regge diagrams for diffractive production in 7*p-interactions. 
Fig. 11 The ratios 1 \ . Experimental points are from ref. p7 



Fig. 12 The ratios F2/AF2N computed from eq. 
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Fig 8 
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